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Abstract 



The rare decay B — > (j)Ks is a well-known probe of physics beyond the Standard Model because 



it arises only through loop effects yet has the same time-dependent CP asymmetry as S — > '4'Ks- 
• Motivated by recent data suggesting new physics in i? — > (pKs, we look to supersymmetry for 

possible explanations, including contributions mediated by gluino loops and by Higgs bosons. 
^ , Chirality-preserving LL and RR gluino contributions are generically small, unless gluinos and 

I squarks masses are close to the current lower bounds. Higgs contributions are also too small to 

explain a large asymmetry if we impose the current upper limit on B{Bs fM'^fj,^). On the 
other hand, chirality- flipping LR and RL gluino contributions can provide sizable effects and while 
remaining consistent with related results in B ^ ipKsi ^Ms, B — > Xg^ and other processes. 
We discuss how the LR and RL insertions can be distinguished using other observables, and we 
provide a string-based model and other estimates to show that the needed sizes of mass insertions 
are reasonable. 
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I. INTRODUCTION 



B (pKs is a powerful testing ground for new physics. Because it does not occur at tree 
level in the standard model (SM), but only via loop contributions, this decay is very sensitive 
to possible new physics contributions to the quark level process h sss. Within the SM, it 
is dominated by the QCD (and electroweak) penguin diagrams with a top quark in the loop 
(see Fig. Therefore the time dependent CP asymmetries are essentially due to B^-B^ 
mixing and thus are the same as those in S — >• ip^s- sm2P^K — sin2/3^x + O(A^) 

Recently the B — > (pKs decay was observed by both BaBar and Belle at the branching 
ratio of ~ (8 - 9) x 10"^ (see Table ^. This is in accord with theoretical predictions 
based on various factorization approximations, albeit there are considerable uncertainties in 
theoretical predictions. Being a penguin dominated process, the theoretical prediction for 
the branching ratio of B ^ <pKs is very sensitive to the so-called effective number of colors 
A'^^, which has been widely used in the old factorization approximations 0. VaryingiV^''' 
from 2 to oo, the predicted branching ratio lies in the range 13 x 10~^ and 0.4 x 10~^ [3|. If 
one uses the improved QCD factorization method developed by Beneke, Buchalla, Neubert 
and Sachrajda (BBNS) 3 (where Nf = = 3), one obtains B{B (pK^) ^ 5 x 10-^ 
somewhat lower than the current world average despite large experimental uncertainties. 
Another approach based on perturbative QCD yields the branching ratio of this decay at 
the level of 10 x 10~® within the SM On the other hand, the CP asymmetries in 
B — > (j)Ks are less model dependent, since they are ratios of branching ratios for B^ and B^ 
decays. Large theoretical uncertainties mostly cancel out between the numerator and the 
denominator (see, e.g., Ref. 0). 

In general, the time dependent CP asymmetry in B^{B^) — > (pKs (or for any common 
CP eigenstates into which both B^ and B^ can decay) measures two independent pieces of 
information 01 : 

, ... ^ ^(B^sit) - <PKs) - r(i?pV(t) ^ <pKs) 

'""^ ^ " r(i?°,,,(t) - <pKs) + r(BO,^,(t) 4>Ks) 

= -C^KCOs{AMt) + S^Ksm{AMt), (1) 



where C^px and S'^x are given by 

_ 1 - \X<pk\'^ AC ^ ImA^i^ , . 

^^'^ ~ 1 I l\ — 12' ^^'^ ^<t>'< ~ 1 , i\ — 12' y"^) 

J- + 1^0x1 t + \A^k\ 



with 



' A{B^<pKs)- 

The angles (3 and 6^ represent the SM and any new physics contributions to the B^-B^ 
mixing angle; we will assume the latter to be small. A{B — ^ 4>Ks) is the amplitude for the 
nonleptonic B decay of interest. If there are several independent channels relevant to the 
same final states with weak phases <^m and strong phases Sm (with m = 1,2,... labeling 
different channels) , the decay amplitude A for a particle is given by 

A = ^a^e+^^'"e+*^'", (4) 



2 



TABLE I: CP averaged branching ratios, S(j,K and C^x in i? ^ (f'Ks from BaBar 0] and 
Belle 3,0 and the SM predictions in the BBNS approach. 



Observable 


BaBar 


Belle 


Average 


SM prediction 


Br (in 10"^) 


8.1 ±0.8 


8.7t|8±1.5 




~ 5 (see text) 


S^Ks 


0.45 ±0.43 ±0.07 


-0.96 ± 0.50 


-0.13 ±0.33 


0.736 ± 0.049 


C(f)Ks 


-0.80 ±0.38 ±0.12 


0.56 ±0.41 ±0.16 


-0.19 ±0.30 


-0.008 



whereas the decay amplitude A for the antiparticle is given by 

A = ^a„e-^^'"e+*'^'". (5) 

m 

Here > is the modulus of the amplitude of the m-th channel. Within the SM, the 
amplitude is real to a good accuracy so that A^/^ = — e~^*^ = X^k- This remains true if 
any new physics contributions are real relative to the SM decay amplitude for B 4>Ks- 
Therefore one needs a new CP violating phase(s) ipm if there exists any significant deviation 
of SipK from S^K- For example, if there is a new physics contribution to the decay with 
amplitude whose modulus is onp, weak phase is ip and strong phase is 6, then 

where r = flNp/osM- In general, the strong phase 6 will affect both S^pK and C^x- However 
if the new physics contributions dominate over the SM ones, the strong phase effectively 
drops out, an approach that is commonly (and mistakenly) used in the literature even when 
r is not large. We will comment further on the nature and source of the strong phases in 
Section lim 

As described before, S'^x ^ sin 2^^^ = S^k within the SM: 

S^K = S^K = 0.736 ± 0.049, and C^k = -0.008. (7) 

That is, the time dependent CP asymmetry in i? — > (pKs should be essentially the same as 
that in B —>■ ipKs- However the BaBar and the Belle collaborations both report a deviation 
from the SM prediction for S^k- As summarized in Table HJ the Belle value for S^k is 3.3(j 
away from the SM prediction, while the Babar value is within la of predictions. (Previously, 
both collaborations had reported values inconsistent with the SM prediction by about 2.7a. 
The average of the two values remains 2.1a away from the SM.) This result, which may be 
an indication of new physics contributions to B — » c^- ft",?, has generated a wave of activity, 
much of it on SUSY contributions io B (pKg H El Q 13, El, H H H 0. The 
direct CP asymmetry in B (pKs is also reported by both collaborations (see Table H]). 
However, their two values have large errors. Measurements of both the direct and indirect 
CP violation are in disagreement between the two experiments, so no firm conclusions can 
be drawn at present. 

However in this paper, we wish to entertain the possibility that there is indeed a deviation 
of sm2j3ii,K from sm2P^K, and that it originates from supersymmetry (SUSY) effects. More 
specifically, we consider the B — > (pKs decay within several classes of general SUSY models 
with i?-parity conservation. We will study two interesting classes of modifications tob —>■ sss 
within SUSY models: 
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• Gluino-mediated b sqq with q = u,d, s, c, b: Such operators are induced by flavor 
mixings in the down-squark sector. We consider all possible combinations of LL, LR, 
RL and RR mixings in the b-'s squark sector. Such contributions do not distinguish 
among flavors of light quarks due to the flavor independence of SUSY QCD. Therefore 
other decays such as B ^ Kn,r]^'^K could be affected as well. 

• Higgs- mediated b sss in the large tan/5 limit (oc tan^/5 at the amplitude level): 
This mechanism is important only for b sss, and not for sec, suu or sdd transitions, 
which makes it an attractive possibility. It would affect (pK and rj^'^K modes but not 
Kt[ modes. 

In this article, we carefully analyze the effects of these two mechanisms on B ^ 4>Ks and 
related observables. More specifically, we consider 

• B ^ 4>Ks '■ its branching ratio and the time dependent CP asymmetries, S^k and 

• Correlations with B — > X^'-f, both its branching fraction and its direct CP asymmetry; 

• Correlation of the Higgs-mediated b — > sss transition with B^ fJ''^f^~ which has 
been, and is being, searched for at the Tevatron; 

• Correlations with B^-B^ mixing coming from SUSY contributions to AMb^ , and with 
the dilepton charge asymmetries An and the time-dependent CP asymmetry in Bg — 
ipcj) (which is proportional to the phase of the mixing). 

The new phase in the b — > sss will affect other CP-violating observables in calculable ways, 
and our explanations can be tested by measuring other quantities as we will discuss. 

Before proceeding, we comment on the following concern: why should there be a large 
deviation in i? — > (pKs but not in related decay modes such as -B rj' Ks,t^K,K^K~ Kgl 
Unlike B (pKs, these decays have SM contributions at tree level while the SUSY contri- 
butions are loop-suppressed. Therefore it is reasonable that only B —>■ (pKs, which is already 
one loop suppressed in the SM, is modified by a significant amount. Because considerable 
hadronic physics is involved, we will not make any more precise statement here, but we do 
not think any contradiction of our analysis of (pKs is implied by the data. The unexpectedly 
large branching ratio for rj'Ks is also consistent with this view. 

We also note that for B —>■ D*^D~* there also is a 2.7a deviation between the SM 



prediction and the data [21|. This decay is dominated by the tree level b — dec transition, 
but its amplitude is suppressed by a factor of A = sin 9^ ~ 0.22 relative to that oi B ^ J / ipK. 
Therefore new physics contibutions at one loop level might have a chance to compete with 
the SM contribution to Bd — > D*^D~* . Possible SUSY contributions involve b — > dqq (with 
q = u,d,s,e) so that the relevant mass insertion parameter is ((51b)i3- This parameter 
is independent of (^1^)23 which affects B (j)K. One could perform a similar study as 
presented below for the Bd — ^ D*^D~*, but we do not pursue that here. 

In the following, we will deduce that LL (and RR) insertions in gluino penguins generally 
provide contributions too small to cause an observable deviation between S^k and S^Ki 
unless gluino and squark masses are close to the current lower bounds. We will also find 
that Higgs-mediated b — » sss is not sufficient to explain the data once we impose the existing 
CDF limit on Bd^s /^/^- However, we find that the down-sector LR and RL insertions in 
gluino penguins can in fact explain a sizable deviation in a way consistent with all other 
data. 
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II. AB = 1 EFFECTIVE HAMILTONIAN 



All approaches to exclusive B M1M2 decays (where Mi and M2 are light mesons) 
use factorization methods at various levels of sophistication. One starts from the AB = 1 
effective Hamiltonian at the renormalization scale fi ~ rrih, which can be obtained from 
the underlying ultraviolet physics by integrating out heavy particles, with the effects of 
hard gluon taken into account by the renormalization-group-improved perturbation theory 
(RG-improved PQCD). 

The effective Hamiltonian for B (pK in the SM can be written as 



cfT 



' p=u,c 



10 



89 



i=3 



+ H.C., 



where \p = V*gVpb with p = u,c are CKM factors, and A„ + Ac + A( = due to the unitarity 
of the CKM matrix. The operators Oi are all those relevant to AS* = 1 hadronic decays of 



a b quark: 










0? = 


{pb)v^ 


-Aisp)v-A 


Of = 


{Pabf3)v-A{SfiPa)v-A 


03 = 


{Sb)v- 


-AY.qi.m)v-A 


04 = 


{sj)p)v-A Y.q{qpqa)v-A 


0,= 


isb)v- 


-A^q{qq)v+A 


Oe = 


{sJ)p)v-A Y.q{W<y)v+A 


07 = 


{sb)v- 


-AT.qleq{qq)v+A 


Os = 


{Sabli)v-A J2q l^qiWajv+A 


09 = 


{sb)v- 


-Aj2qh'l(^l)v-A 


Oio = 


-- {Sabf3)v-AT,g^eq{qf3qa)v-A 


07^ = 


e 

~ 87r2 




Os,= 





(9) 



where a, (3 are color indices. The operators O12 are charged current operators relevant at 
next-to-leading order; 03,4,5,6 are generated by gluonic penguins at leading order; 07,8,9,10 
are electroweak penguin operators, also generated at leading order; Oy^^sg are the magnetic 
and chromomagnetic transition operators, also generated at leading order. The Wilson 
coefficients Ci contain all the relevant information regarding the short-distance physics (and 
possible new physics effects, if any). For simplicity, we will ignore the electroweak penguin 
operators O7...10 whose contributions to B ^ <pKs are roughly ~ 10% within the SM. The 



expressions for Oj's within the SM can be found in Ref. 

The most difficult task is evaluating the matrix element of the above effective Hamiltonian 
between the initial \B) state and the final (M1M2I state. In this work, we adopt the BBNS 
approach |^ to estimate the hadronic amplitude for B — > (pKs- It is inevitable that our 
results will depend on the factorization scheme chosen. The direct CP asymmetry C^k is 
particularly dependent on the method we use for evaluating the hadronic matrix element, 
more so than S'^x- For example, C^k = for naive factorization without one loop corrections 
to the matrix elements of four-quark operators. Including those corrections, or going to the 
BBNS approach, can lead to very large asymmetries. On the other hand, S^k can be 
large and negative in either scheme. We will discuss some of these unavoidable scheme 
dependencies again at the end of Section IVIEI 



III. GLUINO-MEDIATED FCNCS 



In the MSSM, supersymmetric versions of the SM contributions to -B ^ (pKs exist. For 
example, the W-t loop of Fig. ^a) is accompanied a new W-t loop. The flavor-changing in 
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these diagrams is intrinsically tied to the usual, CKM-induced flavor changing of the SM. 
If that were the only new source of flavor physics, we would say the model is minimally 
flavor violating (see Ref. for a consistent deflnition of minimal flavor violation in two- 
Higgs doublet models such as the MSSM). However such a model will not generate large 
corrections to S^k or C^k and so cannot describe the experimental data. (We will return 
to a more-minimal flavor violation when we discuss the Higgs-mediated contributions in 
Section IVH) 

But the general MSSM is not minimally flavor violating. For a generic MSSM a new 
source of flavor violation is introduced by the squark mass matrices, which usually cannot 
be diagonalized in the same basis as the quark mass matrices. This means gluinos (and 
other gauginos) will have flavor-changing couplings to quarks and squarks, which implies 
FCNCs which are mediated by gluinos and thus have strong interaction strengths. In order 
to analyse the phenomenology of these couplings, it is helpful to rotate the effects so that 
they occur in squark propagators rather than in couplings, and to parametrize them in 
terms of dimensionless parameters. We work in the usual mass insertion approximation 
(MIA) [2^ , where the flavor mixing j i m. the down type squarks associated with qs and 
qA are parametrized by More explicitly, 

in the super CKM basis where the quark mass matrices are diagonalized by and V^, 
and the squark mass matrices are rotated in the same way. Here Mg, M|, and are 
squark mass matrices, and m is the average squark mass. Then, the gluino box/penguin 
diagrams generate the QCD penguin operators, Os . g and 03^..^6 = 0-i,,,,(i{L ^ R). Since the 
decay — s> (j)Ks is dominated by the SM QCD penguin operators (an example is shown in 
Fig. n^a)), the gluino-mediated QCD penguins such as that in Fig. ^b) may be signiflcant 
if the gluinos and squarks are relatively light. Similar studies were carried out by Lunghi 
and Wyler ji^ and Ciuchini and Silvestrini (and in the context of a SUSY GUT by 
Moroi [2^ and Causse [lH). Our approach extends these existing papers in the following 
respects: 

• We use the recent BBNS approach to evaluate the hadronic matrix element for B — »• 
(f)Ks {i.e., a'^j^ in Eq. which is irnportant for calculating direct CP violation C^x, 
unlike other papers (except for Ref. (l3|). 

• We consider the LR, RL, LL and RR insertions, whereas Ciuchini et al. consider 
only the LL insertion, and Moroi and Causse consider mainly the RR insertion. We 
will flnd that the Li?, RL contributions can dominate, while the LL, RR insertions 
are too small to affect B —>■ (pK signiflcantly unless the gluino and squark masses are 
close to the current lower bounds. 

• We also study the Higgs-mediated contribution (Fig.^c)), the amplitude of which can 
be enhanced by tan^ /5 for large tan j3. Naively this contribution would be a natural 
flavor-violating candidate to provide a non-SM enhancement of the ss flnal state. 

Note that the gluino mediated QCD penguin operators are not enhanced in the large 
tan P region, unlike the Higgs-mediated contributions. But the gluino-mediated diagrams 
contribute to 6 ^ sqq for all q = u,d, s, c, b, independent of flavor. Therefore one should 
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s s , s 




t(.^^ ((.^^ h,H,A / ~ 



^ ^ b , hcG^CGGC^ > ^ b ^ — / hcG^CCCG^ > ^ 

W g 3C" 

(a) (b) (c) 

FIG. 1: To illustrate different contributions, we show some Feynman diagrams relevant to 6 ^ sss: 
(a) the SM, (b) the gluino mediated and (c) tan /3— enhanced Higgs mediated contributions. 



check that these new sources of flavor changing do not contribute too much to other charmless 
B decays such as i? ^ tttt, nK, etc. In this paper we will not deal with this problem 
directly. But we note that there are known mechanisms by which this can be achieved. As 
an example, consider a scenario in which both LL and RR operators contribute to 5 — > 4>Ks 
with nearly equal amplitudes. Because vector mesons (0, p, ...) and pseudoscalar mesons 
{n,K,T]^^) have opposite parity, the gluino loop effects appear as (5^^)23 + (^r_r,)23 for the 
VP modes and (5^^)23 — {Srr)23 for the PP modes, respectively. So one can suppress the 
gluino contributions to the PP modes by simply assuming (5^^)23 = {^rr)23, even if LL 
and RR contributions are sizable [l^ . 

There are additional issues that arise in comparing B (pKs to B ^ Tin and others. 
For example, there are multiple diagrams contributing to processes such as -B ^ tttt, unlike 
the case oi B —>■ (pKs- This makes the inclusion of SUSY effects in these processes either 
irrelevant or highly uncertain. We will take an indirect approach to this problem by requiring 
that the total branching ratio of 5 — > 4>Ks is consistent with observation after including 
SUSY effects. By demanding that the SUSY contributions to B ^ <f)Ks are of the same 
order or less than the SM contributions, we can safely assume that the SUSY contributions to 
related processes are also close to experiment. Given the uncertainties inherent in the BBNS 
or any other factorization scheme we feel that this is an appropriate approach. However it 
also means that the size of the strong phase will play an important role in our calculations, 
another reason for using the best technology available today, namely BBNS factorization. 

The relevant Wilson coefficients due to the gluino box/penguin loop diagrams involving 
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the LL and LR insertions are given (at the scale ~ mw) by [2 



^SUSY 



^SUSY 



(jSVSY 



^SUSY 



a: 



r 



SUSY 



^I2s 

Snag 



BAx 



-B^ix) - —PAx) 



—BAx) + —B2{x) - —PAx 



:P2{x\ 



5i(x) + ^-B,{x) + iPi(x) + \P2{x)^ 

1P1(X) + ^P,(X)) (5^,),3 



-Ib,{x) + ^B2{x 



23 



^/2GFm^\t 



{Kl)23 ( \m,{x) - ^M4(x) 



4Pi(x) 



Lijj23 
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(4Pi(a;) -9x"iP2(x)) 



(10) 



where x = {rrig/rh)'^, and the loop functions can be found in Ref. [28[. In the presence of 
the RR and RL insertions, we have additional operators Oi=3.. .6,77,83 that are obtained by 
L ^ i? in the operators in the SM. The associated Wilson coefficients Ci=3.. .6,77,89 are given 
by the similar expressions as above with the replacement L ^ R. The remaining coefficients 
are either dominated by their SM contributions (6*1,2) or are electroweak penguins (C7...10) 
and therefore small. 

A remark is in order concerning the above expressions. These Wilson coefficients differ 
from the 1996 results in Gabbiani et al. [2^ 'va. (i) the signs of the (5^/2)23 contributions in 
Cy^ and Csg, and (ii) the loop function 9x-B2(a;) in G^g should be corrected to 9x^^i?2(a^) as 
given above js^]. These differences are not important for numerical results, however. 

The decay rate for B <pKs is given by 



T{B ^ 0/^5 



327r 



B^K\2 I „cfr I 2 



where 



efl 
''(t>K 



p=u,c 



(as + a^ + 05) - - (ay + 09 + a?o) 



(11) 
(12) 



and A is the magnitude of the 3-momentum relative to m^. In the numerical analysis, we 
use frp = 237 MeV and Ff^^ = 0.38. The a^'s are given in terms of the Wilson coefficients 
Gi times various nonperturbative hadronic parameters: 





= C3 + 




as 


Nc 


p 
a\ 


= C4 + 


G3 

— + 




= C5 + 


G, 


05 





C4 G po^s 
Cs Gpas 



GpOts -^1,2 



N.AnN, ^ ''''' 47r N, 
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FIG. 2: A typical final state interaction diagram which generates a strong phase in i? ^ 4>Ks- At 
long distances, the weak phase is contained in the 4-fermion vertex shown as a small square. At 
short distances, that vertex resolves to SUSY diagrams like those studied in this paper. 

where Cp = {N^ — l)/2Nc and Nc = 3. The various V^'\ H^'^ and P terms represent the 
nonperturbative hadronic parameters at the heart of the BBNS calculation and are described 
in Ref. [j. 

Besides generating the amplitudes necessary for calculating physical rates, the BBNS 
approach also provides us with the appropriate strong phase. We will digress for a moment 
to talk about the importance of this piece in the calculation. 

In order to have nonzero CP asymmetries, we need at least two independent amplitudes 
with different weak and strong phases. In the SUSY models we are considering, the weak 
phases reside in the complex mass insertion parameters, the (5's, and appear in the Wilson 
coefficients Cs^,,,^ and Cgg'- see Fig. ^a)-(c). These weak phases are odd under a CP 
transformation. On the other hand, the CP-even strong phases arise from scatterings among 
the final state particles. At the parton level, they are generated by gluon exchange between 
various quarks participating in the B decay processes. In the BBNS approach there are 
four classes of diagrams that can generate strong phases: vertex corrections, penguins, hard 
scattering with spectators and annihilation diagrams. We show a vertex correction diagram 
in Fig. 121 for illustration, referring the reader to the BBNS papers j3| for the details. The 
strong phases are then encoded in the V^'\ H^''> and P terms that appear in the definitions 
of the tti parameters above. 

Before the BBNS approach, the strong phase was generated by the so-called Bander- 
Soni-Silverman (BSS) mechanism, in which the one loop diagrams involving the 4-quark 
operators develop a strong phase consistent with unitarity. However, there is an ambiguity 
in this approach regarding the momentum flow into the quark loop, and one usually made 
the ad hoc assumption that the average through the virtual gluon is (g^) ^ ml/ 4. This 
ambiguity disappears in the BBNS approach, where one has a definite relation between the 
loop momenta of the quark and antiquark and the meson lightcone wavefunction. In partic- 
ular, the contributions of the b sg operators Ogg and Ogg are obtained in an unambiguous 
way at least in the leading order. 
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IV. B Xs-y 

It has long been claimed that the most stringent bounds on the (5^^)23 insertion comes 
from 6 — > 57 (or i? ^ ^sl)- Because the measured value is close to the SM prediction, there 
is little room for new physics lurking into b ^ s transition. This statement is particularly 
true for LR mixing because the SUSY contributions interfere with the SM ones at the 
amplitude level (in Cy^ and Cgg) and further because the SUSY contributions pick up an 
tTig/rrib enhancement relative to the SM. The LL insertion also contributes to C-jj^sg but 
lacks the gluino mass enhancement. Thus LL insertions can be very large while remaining 
consistent with b 57, while LR insertions are constrained to be small, of O{10~^). 

The contribution of the RL and RR insertions to C^j^sg are suppressed by ms/mb where 
rus is the strange quark mass and so we ignore them. However they make unsuppressed con- 
tributions to Cj^^sg, but because these operators do not appear in the SM, their contributions 
cannot interfere with the SM, nor can they generate CP violation. Thus the constraint from 
B Xg'y on RL and RR operators will always be somewhat weaker than that on LR and 
LL operators respectively. 

We will impose a rather generous bound for B Xs'j, 

2.0 X 10"^ < B{B X,7) < 4.5 x 10"^ 

in order to take into account the theoretical uncertainties in the NLO calculation of the 
SUSY contributions. (We have also checked that our results for S(pK do not change very 
much if we narrow the allowed window for the B — > Xsj branching ratio. It will be relatively 
straightforward to see the effect of narrowing the window in the figures we present in Sec- 
tion |V11) The CP-averaged branching ratio for B Xg'j in the leading log approximation 
is given by 



B{B ^ X,7) 
B{B X^eu) 



cb 



2 



'C,,{mb)\'+\Cr,{mb)\' . (13) 



7rf{z 



where f{z) = 1 — 82 + 82;^ — 2;^ — 122;^ log 2; is the phase space factor for the b c semileptonic 
decays and = 137.036. Neglecting the RG running between the heavy SUSY particles 
and the top quark mass scales, we get the following relations : 

C^^irrib) ~ -0.31 + 0.67 C^^^sy^^^^ ^ g gg ^^susy^^^^^ 

Csgimb) ~ -0.15 + 0.70 CfusY^^^^_ (^4) 

In most of our analysis ol B Xg'y we will assume that the only contributions are those 
of the SM plus the gluino-mediated SUSY loops. It is well known that the charged Higgs 
and chargino sectors in the MSSM may contribute to i? ^ ^^7 with strengths equal to or 
greater than the SM piece. We will ignore this possibility except in Section IVI C 11 where 
we will consider a particularly interesting limit in which all but the gluino loops cancel out 
among themselves. 

A new CP-violating phase in ((^^^)23 will also generate CP violation in B — s> ^^7. In order 
to have a nonvanishing direct CP asymmetry, one needs at least two independent amplitudes 
with different strong (CP-even) and weak (CP-odd) phases. In S — >■ Xg'j, strong phases 
are provided by quark and gluon loop diagrams, whereas weak phases are provided by the 
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CKM angles and (^1^)23- For conventional models with the same operator basis as in the 
SM, the resulting direct CP asymmetry in S — > Xs^ can be written as 



A'cp'"(m %) = TTT^ [1.23 Im - 9.52 Im [C^gC;^] + 0.10 Im {C2CQ] . (15) 

Since SUSY contributions to C2 are negligible, we use C2{mh) ~ C|^(m;,) ~ 1.11. In the 
presence of operators Oi with opposite chirality, Eq. (fT5|l is modified according to Ref. jHf . 
Within the SM, the predicted CP asymmetry is less than ~ 0.5%; any larger asymmetry 
would be a clear indication of new physics. 

What we currently know about A^'p'*'^ comes from CLEO 132]: 

= (-7.9 ± 10.8 ± 2.2)(1.0 ± 0.030)%. (16) 

This is clearly not a very strong constraint at present, but should become significantly more 
important in the near future. 

Despite the importance oi B —>■ X^'-f for constraining (5^^)23, we will find in the following 
sections that the B —>■ (j)K branching ratio itself provides a constraint on the LR and RL 
insertions which is every bit as strong as i? — > ^^7. 



V. B^-B'^ MIXING 

Gluino-mediated box diagrams will also affect B^-B^ mixing, which is phenomenologi- 
cally very important. The most general effective Hamiltonian for B^-B^ mixing {AB = 2) 
and its SUSY contributions have been calculated in Refs. [331 l33 and we do not repeat 
that discussion here. Following the explicit formulation of Ref. j3J|, we calculate AM^ and 
the phase of the B^-B^ mixing, 13s. For the S^.s meson decay constants, we have assumed 
fe^ = 200 ± 30 MeV, and 

fp Bn 

f ^ = 1.16 ±0.05. 

Currently, only the experimental lower limit AMb^ > 14.9 ps~^ is known [3^. Also we 
assume that the SUSY effects on B^-B^ mixing are reasonably small jl^ (z.e., Od in Eq. IHl 
is nearly zero). According to the analysis by Ko et al. ji^], the CKM angle 7 can still be in 
the range between —60° and 60° for the LL insertion (5^2.)i3 7^ without any conflict with 
the measured AM and sin2/3^/^g. Then these effects will only appear in as 2{[3 + 6*^), 
which is flxed to be sm2(3^K = 0.734 ± 0.054 by data. Furthermore, the SM contributions 
to 5 — s> (j)Ks or B'^-B'^ mixing are independent of the CKM angle 7, so that we don't have 
to worry about the SUSY contribution to B^-B^ mixing. 

In our model, the CP violating phase in the mass insertion parameters (^^5)23 contributes 
not only to S^k, as discussed before, but also to the imaginary part of the mixing, 
Mi2(-Bs), which is real within the SM to a very good approximation. Thus there should be 
some correlation between S^fyK and 2/3^ = — Arg [Mi2(-Bs)]. The observable S^^ = sm2Ps 
plays the same role in B^ ip(l> as S^k does in _B — 4>Ks, and thus we may expect large 
deviations in the time-dependent decays of Bg into ipcp. The SM predicts Ps vanishes (to 
experimental accuracy), so any measurement of a non-zero (3s would be a sign of new physics. 



11 



The new phase in B^-B^ mixing will also appear in the dilepton charge asymmetry of B 
decays that can be measured at hadron colliders (see, e.g., [s^): 

_ iV(£+£+) - Nji-i-) _ N{BB) - N{BB) ^ ^ 

= NiFnTmFF) ~ N{BB) + N{BB) " ^^^^ri^/Mi^). (17) 

In the SM, the phases of M12 and are almost the same, and Aff^ ^ lO""^ (10"'^ for Bd). 
But in the presence of SUSY, a much larger dilepton charge asymmetry may be possible. 
Though SUSY is not expected to provide a large correction to (which is dominated by 
large SM contributions), it can strongly affect M12. In such a case, the dilepton charge 
asymmetry could be approximated as 

/ pSM \ 

The possible range of Au in a classes of SUSY models were studied in Refs. ^^^^ 
follow the formulation of Ref. [s^] in our analysis. We will demand that all solutions obey 
the lower bound AMs > 14.9 ps~^, and for those that do, present their prediction for AM<j. 



VI. NUMERICAL ANALYSES OF GLUINO MEDIATION 

Now we come to the heart of our analysis. We would like to consider whether gluino- 
induced flavor changing can explain the B (pKs data while remaining consistent with 
all other known constraints. And if so, we wish to identify other observables that can be 
used to check our interpretation of the data. We will consider each possible mass mixing 
(5^^)23 for AB = LL, LR, RL, RR one at a time. The technique assumes that only one 
insertion dominates the new physics. However, where relevant, we will explain how a second, 
subdominant insertion could change some of the correlating observables. 

In all of the studies that follow, we will scan over the modulus and phase of the flavor- 
changing mass insertions (|Re,Im (5^5)23! < 1), calculating all observables for a fixed gluino 
and (universal) squark mass of 400 GeV. We will then vary the ratio x = rn-fm^ away from 
unity (keeping rrig fixed) and finally, at the end of the section, discuss the effects of varying 
the gluino mass. We will keep only those points which are consistent with the branching ratio 
for B — >■ Xg'y, and which have AM^ > 14.9 ps~^. Both of these constraints were discussed in 
previous sections. We will also require that the branching ratio for B — > (pKs be less than 
1.6 X 10~^, as calculated in the BBNS approach. This value is twice the observed branching 
ratio, but we feel that the uncertainties inherent in the BBNS factorization scheme may be 
large enough to merit this window. 



A. The LL and RR insertions 

We will begin by considering the the LL and RR insertions, though assuming that one 
or the other appears dominantly. Since the physics effects of these two operators are almost 
identical, we will couch our discussion in the language of the LL insertion and specifically 
point out those few places where the RR analysis diverges. 

A non-zero LL mass insertion ((5^2.)23 generates the QCD penguin operators Os^,,,^^ and 
the (chromo)magnetic operators The B and AMg constraints remove two 
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distinct regions in the {Re (5^^)235 I^i (5^/^)23} plane. The constraint on the branching ratio 
of 6 — >■ 57 removes all points with Re (5^^)23 < —0.5. The constraint that AM^ > 14.9 ps~^ 
removes points where |Re (5^^)231 < 0.25 and 0.25 < |Im (5^^)23! < 0.75. See Fig. El^a). The 
resulting AM^ is indicated by the different colors. For x = 1, AM^ can be as large as ~ 100 

ps"^ 

However, nowhere in the plane are we able to get substantial shifts in S^k- In particular, 
for rrig = rh = 400 GeV we find no points with S^k < 0.5. (We will comment shortly about 
the effect of changing the SUSY masses.) The lack of impact on S^px arises in part because 
of substantial cancellations between Cf^f^ and Cf^^^. This result is consistent with that 
obtained by Lunghi and Wyler [2^. If we lower the common squark mass to x = 3, the 
asymmetries become somewhat smaller and AMg can be as large as a few tens of ps~^. On 
the other hand, for a heavier squark of x = 0.5, AM^ can be even larger than 100 ps~^, 
which is beyond the reach of hadron colliders. 

(On the other hand, if we lower the gluino mass down to 250 GeV with x = 1, S^k can 
shift down as low as 0.05, but this is possible only if the gluino and squark masses are close 
their current lower limits of 200 to 250 GeV. See the Section IVlDI for further discussion.) 

Furthermore there is very little effect on C^x and A^'^'^ , the former varying only between 
±0.1, the latter between ±3% (see the figure). However there are one-to-one correlations 
among S^k^ C^k and A(^^^ for an LL insertion, so measuring these correlations would 
provide strong evidence for the model. However the experimental precision achievable in 
the near future will make this a very difficult task. 

The lack of signal in B —>■ (pKs and B Xs'y does not mean that an LL insertion is 
without signature. It is clear from Fig. |2fe) that there can be sizable effects in the dilepton 
charge asymmetry Au, with asymmetries as large as an order of magnitude above the SM. 
But the largest effect is in the Bg sector, specifically B^-B^ mixing and the CP-violating 
phase l3s measured in Bs — »• ip(f). Examining Figs. inta),(e) and (f), it is remarkable that 
SUSY could drive AM^ to values as large as 100 ps~^, and could shift the CP asymmetry to 
any value | sin2/3s| < 1. Thus an LL insertion may be better probed in B^ decays than in 
Bd decays such as -B — > (pKs- 

In summary we conclude that the effects of the LL insertion on i? ^ are insignificant 
unless the squarks and gluino are very close to their current experimental limit. For squarks 
and gluinos with masses of several hundred GeV or more, the LL insertion alone is inadequate 
for explaining the current data on Si^k- However, the LL insertion plays its strongest role 
in the Bs system, where it can contribute sizably to B^-B^ mixing even for heavy squarks 
and gluinos. Thus there remains the possibility that AMg and sin2/3s will be quite large 
while B — > (f)Ks and B X^^ remain close to their SM predictions. 

The physics of the RR insertion is identical in almost every way to that of the LL 
insertion. The prime exception is the b ^ s'j constraint, which is weaker because the RR 
contribution does not interfere with the dominant SM amplitude for B — > ^^7. 

Therefore, the entire region of {Re (5^ij)23, Im (5^/{)23} with 1(5^^)231 < 1 is now allowed 
apart from the small areas in which AMg < 14.9 ps~^: see Fig. IH^a). In Fig. |3J we show 
the plots for the RR insertion case with rrig = m = 400 GeV. It is obvious that the RR 
case is essentially identical to the LL case in Fig. El There are two issues which are not 
obvious from the figures. First, as already indicated, the RR insertion contributes very little 
to B ^ Xg'-f. Thus we have omitted a plot of the prediction of the direct CP violation in 
B Xg'-f, A^(^'^'^ , since it would differ little from the SM. Secondly, for light gluinos and 
squarks, near their experimental limit, it is possible to get slightly lower values of S^i^k for 
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FIG. 4: (a) The allowed region in the (Re ((5^j:j)23, Im ('^^/j)23) plane, and the correlations between 
(b) S(i,K and B{B 4>Ks), (c) S^i^- and C^k, (d) S^i^ and sm2f3s for = 400 GeV and x = 1. 
Different shades are used for different ranges of AMg. 

the RR insertion than it was for the LL; this is because the LL insertions required to further 
push down S^k are inconsistent with B Xg'y. 

A few general remarks on the LL and RR insertions are in order. First, we note that 
one requires a large 2-3 mixing in order to obtain any observable effect in S^k- But for such 
large mixing, our mass insertion approximation is no longer valid. One must then consider 
the full vertex mixings in the di — dj — g couplings with squarks in the mass eigenstates, 
as done in Ref. jl^. In this case, the superGIM suppression is less effective and the loop 
functions are enhanced, which could lead to somewhat larger effects on i? — >■ (pKs- 

If one considers a large 2-3 mixing in the LL or RR sector, an important constraint from 
B^-B^ mixing, which is proportional to (51b)i3, may also arise. Since the SM contribution 
to AMrf is proportional to (VtftV^^)^, the gluino loop contribution with light gluinos/squarks 
can easily dominate the SM, unless the RR or LL insertion for b ^ d transition is very 
small. In Ref. ji^l, two of us discussed the allowed region for the (5^^)13 insertion, imposing 
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the measured AM^, sm2j3^K and B — s> X^'j- It was found that one needs {S'Il)i3 < 5 x 10~^ 
even if we assume that the CKM angle 7 is a completely free parameter. This should be 
compared to (5^^)23 ~ 1 which is needed in order for the LL insertion to explain the observed 
shift in SfpK- It would be difficult, although not impossible, to build a flavor model where 
this is the case. 

B. The LR insertion 

The case of an LR insertion of the form s^j^b^ is very different from that of either LL or 
RR. While LL, RR insertions can be 0{1), the parameters (5^b)23's must be small in order 
to avoid excessively large FCNC amplitudes in charmless nonleptonic B decays. But this 
also implies that even small ((^1^)23 's can lead to observable effects in the B sector. The 
analysis of the LR insertion is particularly simple, in part because it does not contribute to 
the penguin operators C3...10. But its contributions to the (chro mo) magnetic dipole operators 
can be quite large because it breaks the SM chiral symmetry in the 6-quark sector allowing 
dipole operators to occur without the usual suppression, replaced instead with a SUSY- 
breaking rrig insertion. Thus, as previously noted, the LR insertion is much more strongly 
constrained by -B — > ^^7, roug hly |('^^/^)23| < 10"^. However the same chiral symmetry 
breaking also ensures that the LR contributions to Cgg are likewise large, and so there may 
still be a residual effect on B ^ <pKsi both in its branching ratio and CP asymmetries. 
(This last effect was also noticed in Refs. [itI l2ol|). 

In Fig. Efa), the allowed region for {Re(5^^)23, Im(5^j:j)23} is shown for m = 400 GeV 
and X = 1. The dark regions are consistent with the B{B ^ ^s7) and AMg constraints, 
but are not constrained by the bound on the B <f)Ks branching ratio. We impose this 
constraint by hatching the regions in which B{B <f)Ks) > 1.6 x 10"^. Note that there are 
large regions of parameter space which are consistent with 6 ^ 57 but not B{B ^ 4>Ks), 
as we advertised. And though (5^^)23 is constrained to be < 10~^, it is still very important 
for B — >■ 4>Ks, affecting both its branching ratio and the asymmetries S^px and Ci^k by 
significant amounts. In particular, it is encouraging that the branching ratio for B — * (pKs 
can easily be enhanced compared to the SM, moving it closer to the experimental data. 

Fig. El^b) indicates that the branching ratio for B — *■ (pKs can in fact be greatly enhanced 
relative to the SM case. Also note that S^k can take any value between —0.6 and 1, even 
after we impose B{B (pKs) < 16 x 10~^ as a constraint. This is because the SUSY 
amplitude has a size similar to the SM amplitude so that the resulting branching ratio is not 
too different from the SM prediction. But it has very different strong and weak phases from 
the SM amplitude so that S^k and C<^x can vary greatly with respect to SM predictions. 

In Fig. E^c), we show the correlation between S(j,K and C^k- Note that we can get 
negative S^j^x, for which C^k also becomes negative. Conversely, a positive C^px implies 
that S(f,K > 0.6 in our model. Thus the current data on S^k argues strongly in favor 
of a negative C^k if the LR insertion is dominating the new physics contributions. This 
correlation between S^k and C^pK will be tested at B factories in the near future. 

The correlation between S^k and the direct CP asymmetry in B Xs'~f (= A^'^^) is 
shown in Fig. E^d). We find A^f^^"^ becomes positive for a negative S^k, and a negative 
Af^'^^ implies that S^k > 0.6. However the present CLEO bound on A^'^^ fEq. [T^ cannot 
constrain the LR model very much. 

In Fig. El^e), the correlation between S^k and An is shown. Note that A^ ~ 10^^ within 
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FIG. 5: (a) The allowed region in the (Re ((5^^)23, Im ((5^^)23) plane, and the correlations between 
(b) S^K and B{B 4>Ks), (c) S^px and C^pK, (d) S^k and A^"^"'^, (e) S^k and Au, and (f) S^px and 
sin 2/3^ in the LR insertion case for rrig = 400 GeV and x = 1. The dotted boxes show the 1 a current 
expermental data, and the hatched regions represent the region where B{B 4>^s) > 1-6 x 10~^. 
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the SM so that any appreciable amount of Au will be a clear indication of new physics in 
B^-B^ mixing. For our model with the LR insertion, the deviation of An from the SM 
prediction is small after imposing the B — >■ Xg'y and B — s> (pKs branching ratio constraints. 
Finally in Fig. |31^f), we show the correlation between S^k and sin 2/3,, (the latter of which is 
zero within the SM). Note that sin2/3s is nothing but 5"^,^ in the Bg ip(f) decay, and can 
be measured at hadron colliders if AMg is not too large. However, the resulting sin 2/3s for 
an LR insertion is too small to be observed. 

If we lower the common squark mass to x = 3 for a fixed rrig = 400 GeV, the asymmetries 
become somewhat smaller and AMs can be in the range [15.8, 16.2] ps~^. On the other hand, 
for a heavier squark x = 0.5, the asymmetries can be larger, but we still have AMg ~ AMf^. 
Thus the results are not particularly sensitive to the gluino and squark masses. As the 
masses get very large the SUSY contributions slowly decreases, and are too small to provide 
an explanation for a negative S^k- 

Before closing this subsection, let us comment further on the sensitivity of our results 
to the methods used for calculating the hadronic amplitude. It is clear that the direct 
CP asymmetry Cff^x depends very strongly on the method used for evaluating hadronic 
matrix elements. For example, one finds = in naive factorization (without one loop 
corrections to the matrix elements of four-quark operators), whereas it can take on values 
of C(l) in the BBNS approach. The calculation of S^k is also technique-dependent, though 
less so. Different calculation schemes can produce values of S^k which are different by 
factors of a few. For example, S^k can go as negative as —0.6 in the LR insertion case 
when we use the recent BBNS approach, but when the same parameter space is studied 
using the techniques of Lunghi and Wyler 2^|, S^pK only drops as negative as —0.2. We 



consider the BBNS approach to be the most accurate and well-motivated available, which is 
why we use it, but view the differences between these schemes to be an unavoidable source 
of uncertainty at present. 



C. The RL insertion 

The final case to be considered is perhaps the most unusual. An RL insertion is of the 
form s^^L, coupling the RH strange squark to the LH bottom squark. Our intuition from the 
SM would lead us to believe that such insertions would be small compared to LR insertions, 
naturally suppressed by ~ nis/mb- In SUSY models with minimal flavor violation, this is 
indeed the case. However, once one moves beyond minimal flavor violation there is no reason 
for the RL insertion to be particularly suppressed with respect to LR. Furthermore, it has 
a phenomenology that is quite different than the LR insertions because it does not interfere 
with the SM contributions which are themselves of the LR form. (Because of this, we had 
to extend the BBNS approach to the case with opposite chirality 4-quark operators.) A 
quick comparison of the results for the RL insertion in Fig. [Ulto those of the LR in the last 
section should convince the reader of this easily. 

Though the RL insertion generates a wide range in B{B — > (pKs) like the LR case 
(Fig. injb)), the similarities end there. Because the RL insertions do not directly interfere 
with the SM, the allowed region for (5|3)hl is centered around zero (that is, only its square 
is relevant in most observables) as seen in Fig. Efa). And although both cases can easily 
generate large, negative S'^x and negative C^k, only in the LR case is the correlation 
completely clean. As we pointed in the previous section, an LR insertion with negative 
S(i,K implies negative (but not very large) C^k- For the RL case, negative S^k can be 
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FIG. 6: (a) The allowed region in the (Re ((5^^)23, Im ((5^^)23) plane, and the correlations between 
(b) S^K and B{B — > 4>Ks), (c) S^pK and C^x, and (d) iS^x and AMg for rrig = 400 GeV and 
X = 1. Current data at la level is shown by the dotted boxes, and the hatched regions represent 
the region where B{B — > (pKs) > 1.6 x 10~^. 

accompanied by C^k of either sign (Fig. El^c)). However there is a considerable hole in the 
figure around C,^/^ = 0, which may help in disentangling the structure of the insertions once 
more data is to be had. 

Otherwise RL insertions leave very little other evidence. They do not appreciably alter 
AMg, and they do not generate observable ^^^'^ Thus for a pure RL insertion, one 

expects to reproduce the SM in almost all ways except in B <pKs- In fact, of all the cases 
we have considered so far, the RL case best fits the current data, for precisely this reason. 

There is a particularly interesting variation on the RL theme that we now turn our 
attention towards. 
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1. A special case: the RL dominance scenario 



Thus far, in our calculations of -B ^ Xg'y we have made a simplifying assumption that 
can often be incorrect in realistic models. We have assumed that the calculation of -B ^ X^'j 
consists entirely of a SM piece and gluino-induced SUSY pieces. However low-energy models 
of SUSY usually have several other, often large, contributions to -B — Xs'j, most particularly 
from charged Higgs loops (which add to the SM terms), and from chargino loops (which 
subtract from the SM). It is well-known (see, e.g., Ref. [s^l) that in the supersymmetric 
limit, the SM, charged Higgs and chargino contributions cancel each other exactly, predicting 
a rate for B of precisely zero. If SUSY is broken but the superpartners are relatively 

light, much of that cancellation can still occur. In such a case, the observed rate for 6 — >■ S7 
would be due almost entirely to the gluino loops. 

Ref. js^ has argued for just such an unconventional interpretation of 5 ^ -^s7- Here 
the Standard Model and usual SUSY contributions to Cj^ approximately cancel: C7^(m(,) ^ 
Csg{mh) ^ 0. Then the dominant contribution to 6 — 57 is from the gluino-fe^-S/j penguin in 

C7, and is proportional to (5^x)23- We will call this the i?L-dominance scenario in 6 — > 57. 

Now the constraint from B X^'j plays a completely different role in constraining our 
parameter space. Rather than limiting the size of the RL insertion to be close to zero, it 
actually demands a finite non-zero value for the insertion in order to reproduce the observer 
branching ratio. Thus, it has the effect of constraining the complex parameter (5^l)23 to lie 
in an annulus, just as the constraint from the B — > (pKs branching fraction does (though 
of different center and radius). The resulting allowed parameter space is shown in Fig. [7|^a) 
for rrig = rh = 400 GeV (x = 1). In Fig. Ulp) and (c), we show the correlation of S,j)K with 
B{B — > (j)Ks) and C^k, respectively. Note that S'^x can take almost any value between — 1 
and +1 without conflict with the observed branching ratio for B — > (pKs. Also C^k can be 
positive for a negative S^k, unlike the LR case. In particular, \C^k\ > 0.2 for S^ftK < 0. 
This could be a useful probe for identifying this scenario. Since there is only one diagram 
contributing to 6 ^ 57 in this scenario, one has ^4^"^'*''' = trivially. 

However, as discussed in Ref. js^, if (5^/j)23 is also nonzero then A^'*''^ arises as an 
interference between Cj^ and Csg- Then the RR insertion can contribute more significantly 
to AMg] we will not study this more complicated scenario here because the observables 
will depend crucially on the relative size of the various insertions. For example, the RL 
insertion could be induced by a large LL or RR insertion coupled with an extra SU(2)- 
violating mass insertion. In such a case, there could be large deviations both in S^k and 
B^-B^ mixing. Nonetheless, the correlations between S^k and C^k that we demonstrate in 
this paper depend primarily on the LR or RL insertions, and thus are good probes of new 
physics contributions to -B — > 4>Ks. 

Finally in the allowed region of {Sfii)23, we find AMg < 16.5 ps~^, whereas the SM 
prediction is AMg = 16.2 ps~^. If the pure i?L-dominance scenario is realized in nature, 
B^-B^ mixing will be observed shortly at the Tevatron. On the other hand, both Agi and 
sin2/3s are too small to be observed. 

If we lower the common squark mass to x = 3, the asymmetries become somewhat smaller 
and AMs falls in the range [15.7, 16.7] ps~^. On the other hand, for a heavier squark x = 0.5, 
the asymmetries can be larger, but AMg ^ AMf^. 

In summary, LR and RL insertions of a size ~ 10~^ - 10~^ can induce large shifts in S^k 
and C^K (and ^c"?'*'*^ ^'^^ insertion or for RL + RR insertion) but will not affect B^-B^ 
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FIG. 7: (a) The allowed region in the (Re ((5^£^) 23, Im (5^^)23) plane in the special case of sec- 
tion I VI (TTl and the correlations between (b) S^k and B{B — > (pKs), (c) and C^X) and (d) 
and AM^ for = 400 GeV and x = 1. Current data at la level is shown by the dotted 
boxes, and the hatched regions represent the region where B{B 4'Ks) > 1.6 x 10~^. 

mixing, whereas the LL and RR insertions with a size > 0.2 can induce large changes in B^- 
mixing (both in its modulus and phase) but are still too small in generic models to affect 
S(i)K, C(t>K and ^c'p'*'^ significantly. By understanding how each mass insertion individually 
affects the various observables, it is rather straightforward to go from here to the case in 
which more than one insertion is present. We will not complete that study here. However 
there has been recent work that benefits from multiple insertions that arise in well-motivated 
ultraviolet physics scenarios and we refer the reader there for further details. 
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D. Gluino mass dependence of 



In order to generalize our results and present a more complete picture of SUSY contri- 
butions in B ^ (pKs, it is necessary to study how S^k changes as a function of the SUSY 
mass scale. We have already considered how changes in the squark mass scale, with the 
gluino scale held constant, will affect the observables. Now we turn to the gluino mass. 

In particular, we have calculated S^k as a function of the gluino mass, for x = 1, in each 
of the four insertion scenarios: LL, RR, LR and RL. We have shown the results of our 
calculations in Fig. |H1 In each of the four figures we have shaded in the allowed region of 
S^K as a function of rrig. The different shades represent different predictions for AMg as 
described in the figure legends. It is clear that S(j,K can have a large negative value very 
easily for the LR or RL insertion cases over a wide range of gluino and squark masses. 
On the other hand, in the LL or RR case S^k can have a large negative value only for a 
very light gluino (close to the current bound), and quickly reduces to the SM prediction 
as the gluino mass increases above ~ 300 GeV. For heavier squarks x = 0.5, the window 
for a negative S^k increases slightly, but again the effect goes away quickly as nig becomes 
larger than 300 GeV. Furthermore, this always accompanies a very large AMg (greater than 
100 ps~^). Thus Run II can rule out an RR/ LL explanation for S(I,k if gluinos are not found 
and/or AM^ < 100 ps~^ is measured. 

These results speak directly to the analyses found in Refs. jisl . 16l |. Both of these analyses 



followed the same general lines as ours, but found that RR insertions could sizably alter 
S^K- The reasoning behind these claims is now clear. RR insertions (and also LL to a lesser 
extent) can indeed produce S.^^ < 0, but only for very light sparticles. Specifically, we find 
that TTig must fall below 300 GeV in order to obtain S^k < with an RR or LL insertion. 
Like us, the authors of Ref. find that large changes in S^pK must be accompanied by 
large shifts in AMg, making it unobservable at the Tevatron. Ref. disagrees with this 
finding, though the reasons probably have something to do with each groups' estimations of 
the uncertainties in the BBNS factorization procedure. (See Section lVl El for more discussion 
of this point.) 

One note of explanation is needed in order to understand the graphs, and in particular 
why the LR and RL insertions seem to show no sign of decoupling, while decoupling is 
readily apparent for the RR and LL insertions. The reasoning quite simple. For the RL 
and LR insertions, the leading constraint comes from B — > Xg'y, which decouples at the 
same rate as the new physics contributions to Sf/^K- Thus, as rrig increases, the strength of 
the B — ^ Xg'j constraint decreases, but the ratio of {Sfj^ j^j^) 23/ m-g remains constant. On the 
other hand, there are no strong constraints on the LL and RR insertions, so even for light 
gluinos we have allowed values of {S'[l,rr)'23 as large as 1; as the gluino mass increases, we 
are not free to increase the insertion any further to offset the mass suppression. Thus the 
effects on S^k fall as l/m~. 



E. Hadronic uncertainties 

Finally, we should address the hadronic uncertainties within the BBNS approach to fac- 
torization. These uncertainties appear to be the origin of the qualitative differences that are 
found between our results in the previous subsections and those in Ref. [l6| . 

BBNS factorization fails two possibly important ways: it does not properly account for the 
subleading power corrections coming from annihilation diagrams, nor for those coming from 
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FIG. 8: S^K as a function of rrig for x = \. The shading represents AM^ in the ranges [14.9,20], 
[20,40], [40,60], [60,80], [80,100] and [100,cx)] ps"\ respectively. 



hard scattering with a spectator quark. Each of these contributions should be subleading, 
but each involves an integral that blows up in the deep infrared (z.e., Jg dy/y). Thus these 
"subleading" corrections are effectively infinite. In order to control the size of the integral, it 
is necessary to cut it off. We follow the original work of BBNS in parametrizing the integrals 
as A = (1 + pe*'^) \og{mB / /^h) where \h = 500 MeV. For the results we have presented in 
this paper, we have taken A = 0; that is, we have assumed that the subleading corrections 
are very small. 

However, we now wish to see how our results depend on the size and phase of Ahard and 
Aannih- Again followiug BBNS, we vary p and (f in the ranges < p < 1 and < (p < 27r. 
(Each Ahard.annih has its owu p, (f which are varied independently.) In Fig. IHlwe consider the 
RR insertion, since most recent works have emphasized this. For this case, the most negative 
value of S^K occurs at (523)-RR = ^ ~ h for fixed nig = 250 GeV and x = 1. With this value 
of {S23)rr, we plot 5'^;^ as a function of nig. In Fig. Efa), who show the result without the 
subleading corrections. In (b), we turn on the hard scattering corrections only; in (c) we 
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do likewise for the annihilations corrections. Finally in (d) we allow both Ahard and Aannm 
to vary over their entire "allowed" ranges. Note that the annihilation diagrams generate 
much larger hadronic uncertainties than do the hard scattering amplitudes. Theoretical 
uncertainties decrease quickly for heavier gluinos, and there is no possibility to get a negative 
S^K for rrig = rh = 400 GeV within the mass insertion approximation. 

In the work of Ref. the authors did claim to find negative S^pK without generating 
a large B^-B^ mass difference. In order to avoid a large AMg, a smaller value for the 
RR insertion was used. Yet our results would indicate that a small RR insertion will not 
generate highly negative S^k- The solutions to this paradox appears to lie in how the authors 
of Ref. 0] treated the hadronic uncertainties in the BBNS prescription. In particular, they 
allowed phard and pannih to take on values as large as 8. It is not surprising that this 8-fold 
increase in the uncertainties allows for a much wider range for S^px. On the other hand, at 
these extremes, the "subleading" corrections are no longer subleading, but are every bit as 
large as the leading terms in the BBNS scheme. For our analysis we have felt it important 
to keep our theoretical uncertainties at the level of those proposed by BBNS. It seems to us 
that doing otherwise it tantamount to claiming that the BBNS prescription is not a good 
one. We have taken the opposite view in this paper. 



VII. HIGGS-MEDIATED FCNCS 

At large tan /3, FCNCs can also be mediated by exchange of neutral Higgs bosons, as found 
in 5 — ^ Xsl^l~ js^, 5s, d — ^ ji^l and r — >• 3/i j4l|, etc. In this paper we consider the 

same mechanism leading to 6 — sss. Unlike the gluino mediated amplitudes, which generate 
flavor independent amplitudes for h —>■ sqq, the Higgs exchange diagrams generate flavor- 
dependent amplitudes since the Higgs coupling is proportional to the Yukawa couplings. 
Therefore b —>■ sss can be modified by a significant amount, whereas b — suu or 6 ^ sdd 
remain essentially unchanged. Thus decays such as B ^ (pKs,'!]^''^ K can be affected by 
Higgs exchange, whereas B — > nK, nn modes are not. 

In the presence of Higgs-mediated FCNCs (Fig.[TJc)), new operators O^R^Lj^s and 0^r^l}^p 
appear in the AB = 1 effective Hamiltonian: 

^ ifeff - ^KVtb [Cr,sOr,s + Cr,pOr,p + {R^ L)] (19) 

where 

Or^s = (sLbp) (ss), Or^p = (sLbp) (s75s) (20) 

where the first index R ot L denotes the chirality of the initial b quark in the B meson, 
and color is conserved pairwise. In the MSSM with minimal flavor violation, the Wilson 
coefficients Cl^s and C^^p are suppressed relative to Cr^s and Cr^p by m^/mfe, and will be 
ignored in the following. The explicit form of Cr^s{p) due to the neutral Higgs exchange is 
given by 

a f tan^ p \ f mbm^mtiJ,\ sin26'f 2 2 2^ /oi\ 
Crs — —Crp ~ — 7j , ^„ , ^„ f (nif , m,- , a ) . (21) 

where the loop function f{a,b,c) can be found in Refs. ji^], for example. Note that the 
Or^s — Or^p operator is equivalent to — after Fierz transformation (but without summing 
over all flavors, only with ss). 
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(c) Ahard = 0, Pannih < 1 (d) Phard < 1, Pannih < 1 

FIG. 9: Range of predicted in the BBNS approach for several parametrizations of the sub- 
leading uncertainties, as a function of nig. In (a), all subleading uncertainties are taken to be zero. 
In (b) and (c), hard scattering and annihilation corrections are each included. In (d), both are. 
Plots take ((Jfa)/?/? to generate most negative S^k, for rrig = 250 GeV and x = 1. 

In a model with minimal flavor violation (as defined in Ref. [i^) the phases of the SUSY 
contributions would match those of the SM and so no new source of CP violation would 
be generated. However, if we extend minimal flavor violation to allow arbitrary phases (a 
definition used by other authors), we can allow our Cr^s and Cr^p operators to be complex 
and contribute to S^k. Therefore we will assume that the new operators have complex 
coefficients with 0(1) phases. 

However the Wilson coefficients Cr^s and Cr^p also generate Bg — > fi'^fi^, which has 
been and is being searched for actively at the Tevatron. The upper limit on this decay from 
CDF ji^] during the Tevatron's Run I is 

B{Bs /i+/i~) < 2.6 X 10^^ (22) 
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The relevant effective Hamiltonian is 

H^siB, ^ = Vt,VZ [CnA'^tb^R) if^f^) + CnAslb^R) (/^75/i) + (R ^ L)] . (23) 

Note that the Wilson coefficients Cfj^s{P) here are essentially the same as those that appear 
in the effective Hamiltonian for b — ^ sss, up to a small difference between the muon and 
the strange quark masses (or their Yukawa couplings). We find that the current upper limit 
fl22|) puts rather a strong constraint on Cr^s{p)'- 

\Cr,sI\Cr,p\< 0.005, 

since 

B{Bs ^ fifi) oc \Cn,sf + \CB.,pf • (24) 

We then find that S^k cannot be smaller than 0.71 for such small values of Cji^s{P) within 
the naive factorization approach. Thus it appears impossible to explain the large deviation 
in S(f,K with Higgs-mediated b — >• sss alone. 

Recently the authors of Ref. have re-examined this issue in models with non-minimal 
flavor violation. Among other things, they allow a large (6LSij)(ss) operator, though one 
would normally expect such an operator to be suppressed by iris/ nib with respect to the 
{bRSL){ss) operator. The authors of Ref. jl^ argue that one can explain the data on S^jjk 
without violating the CDF bound on Bs if one allows the original operator and its 

parity conjugate to both be large and roughly the same size. 

Since our explanation for S^k involves non-minimal flavor violation, we should in principle 
allow similar violation in the Higgs operators. In the general case the branching fraction for 
Eg fifi depends on 

B{Bs ^ ^/i) oc \Cr,s - Cl,s\' + \Cr,p - CL,pf (25) 

rather than Eq. ()24|) above. It would then appear that our limit applies to this difference 
rather than the coefficients individually [i^. The authors of Ref. jl^l make use of this to 
find a region of parameter space consistent with negative S^px- 

However, this is not the complete story. Because the couplings and masses in the MSSM 
Higgs sector are highly constrained {e.g., rriH — rriA), the scalar and pseudoscalar coefficients 
are tied to one another. In particular, one finds that 

Cr,s - -Cr^p 
Cl,s ^ +Cl,p 

where the approximate equality becomes exact for rriA ^ friz- Thus Eq. ()25p can be rewrit- 
ten as: 

B{Bs ^ M oc \Cr^s - CL,sf + \Cr,s + CL,sf = 2 {\CrM' + 1^^,51') . (26) 

Our argument from the CDF bound applies as before, and we see no way to use Higgs 
mediation to generate S^k < 0, at least in this way. 
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VIII. MOTIVATING {5iRRL)2'i OF 0(10-2) 

In this section, we want to provide some motivation for (5^^)23 ~ IQ-^ that we find 
phenomenologically is needed to generate S^k- Conversely, considerations such as those in 
this section show how one might use data to point toward classes of string-based theories. 
At this stage, we are not advocating any model of {S'[j^)23 ~ IQ-^, but rather showing that 
(5^^)23 of the needed size can be theoretically plausible. 

Before going to a specific construction, we write down the generic expression for the 
flavor parameters (6) lr in a supergravity Lagrangian (which is the low energy effective field 
theory derived from a string theory model). Under the assumption that either the Yukawa 
coupling does not depend on the moduli fields, or it is a function of moduli times some 
constant proportional to the Yukawa coupling itself, we can write the trilinear couplings 



as 



ii 



4 cx (Al ■ Yf + Y-f ■ An),^, (27) 

where Al and Af> are 3x3 diagonal matrices, and f = u,d. We also neglect a term which 
is proportional to a constant times the Yukawa matrix which will not affect the prediction 
of {S)ij. The mass insertions are defined as the ratio between the off diagonal elements 
and the diagonal ones of the squark mass matrices in the superCKM basis. The superCKM 
basis is defined by rotating the quark superfields by V/ and V^, where V/F-^V^^ = Idiag- 
Then we obtain the expression for the relevant LR mass insertions in the general MSSM as 

ihRh,^, oc mjiA^ - A^)iVLUVlh2 + mf^{A^ - A^)iVLUVlh 
+ m{ {A^ - Af){VnMV^),2 + m{ {Af - A^){VnUV^U 

(28) 

where Af'^ are model dependent parameters of order the gravitino mass, 171^/2. The most 
important feature of this expression is that the sizes of various terms are proportional to SM 
fermion masses m{ ■. Some of the immediate conclusions one can draw from this observation 
are 

• (5lr)i2 < (5Li?)l3, (5Li?)23- 

• 5^/6^ ~ m"/m^. In particular, {5'^R)2'i/ {^tiil'i'i ~ rrit/mb. 

Note that these are relations obtained at some high energy scale. However, the dominant 
RGE runnings of the trilinears are diagonal in the superCKM basis. Therefore, those relation 
also hold approximately at low energy scale. 

We now turn to a possible string theory motivated scenario where those features can 
arise. For definiteness, we focus on the models constructed from a intersecting D5-brane 
setup. Suppose we have two stacks of D5-branes intersecting at 90°. There are six compact 
dimensions which are grouped into three complex pairs labeled by i = 1,2,3. D-branes 
are objects on which open strings can end. Open string sectors in this type of scenario are 
classified by (i) which D-brane it ends on, and (ii) which complex dimension it moves in. 
We will have the following open string states ji^ 

• Cj ': open strings which end on the 5^ brane and move along the jth complex compact 
dimension. 

• C*^'^^ : open strings which start from the 5j brane and end on the 5,- brane. 
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The string theory selection rules permit two type of large couplings (provided they are gauge 
invariant): Cf'C|'Cf' and Cj'C^'^>'C^'^>' , where %^ j ^k. 

The next question is how to embed the MSSM matter content into the this scenario. 
Since we need some large mixing between the last two generations, one obvious choice is 
to embed them identically, z.e., assign them to the same open string sector. However, this 
type of construction has a problem because it also predicts = A^ which is not desirable. 
Therefore, we have to embed the last two generations into different open string sectors and 
yet give them a large mixing. 

Suppose we have the following embedding of the MSSM matter content into open string 
sectors, adopting the notation that Fi and Fi are the ith generation left and right handed 
quarks, respectively. 

• The Higgs field: h = ^(Cf + Cl'). 

• The quarks: 





/^5i 

— ^2 ' 


F3 




F2 


= c!\ 


F2 




Fi = 




Fi 





From the string theory selection rules, we can then determine the Yukawa coupling matrix 
to be 

(30) 




We can also work out the soft supersymmetry breaking terms |4y]. In particular, the 
trilinear couplings are 

A 1 ( + -^1 + 

- -— (Xo + X3)-r+ X1 + X3 




1 

1 



+ y\ 1 h (31) 



where Xj's are complex parameters satisfying Yl^=o l^iP = 1- Notice now we have written 
the trilinear couplings in the form of Eq. ()27|). We have singled out the part which is the 
Yukawa matrix multiplied by a diagonal matrix proportional to the identity since this term 
is diagonal in the super-CKM basis and does not contribute to flavor violating processes. 
We also have the following expressions for soft masses 
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(32) 
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Using Eq. ()28p . our ansatz for the Yukawa textures and the definition of the mass insertion 



parameter, we then have the following estimate for the size of {Sfj^ 



1 23- 



Clearly, this can be as large as of order ~ 0.01. It can also have a quite nontrivial phase 
structure since Xq and are in general complex parameters. 

Notice however in this particular implementation, we will generically have (5^/2)23 ~ 
since m|,^ = m'p^. We can also choose (5^2^)23 to be close to zero by setting \Xo\ ^ IX^]. 

As remarked above, this model suggests (5^/2)23 ~ ^ {^lr)23, so (5^^)23 ~ 0.008 suggests 
{61^)23 ~ 0.3. Such a large value is not in confiict with any data on i? — > Xg'j or 5 — >■ 
XJ~^l~ , which give essentially no constraint on (52_r)23 Also (52_r)23 ~ 0.3 can generate 
a fiavor changing top decay t eg only up to (6.0,8.5,9.2) x 10~^ for x = (0.5, 1,3) and 
rrig = 400 GeV, which are still below the threshold for observability (~ 10^^) at future 
colliders such as LHC and NLC Another prediction of this model is (5lr)i2 ~ 5 x 10~^, 
which could dominate e /ex through the chromomagnetic s ^ dg transition. 

Let us also comment on the size of (5^/^)23 in our current embedding. From our formula 
for the trilinears, we see that (5^^)23 = ('^lr)32 is determined by (A^g — )VrV^. In our 
current embedding, we have A2 = Af^ since the last two generations of the right-handed 
(s) quarks have the same embedding. We can then use the unitarity of Vr to simplify the 
expression of (5rj^)23 to (5^^)23 oc ^^y^(VR)3i(VR)2i which is suppressed by higher powers of 

A. However, we can realize a larger (5^/^)23 by simply switching the open string embedding 
of the last two generations of left-handed and right-handed (s) quarks, ^2^3 and ^2,3. Then 
we can have a sizable (5^^)23 ~ 0.01 while (5^^)23 is suppressed by the higher power of A. 

Before closing this section, let us comment on the mass mixings/insertions that would be 
generated by the RG running from the high energy scale to electroweak scale. If we assume a 
universal boundary condition for scalar masses at the reduced Planck scale M^, ~ 2.4 x 10^^ 
GeV, and all {SAB)ij are zero at the high energy scale, nonnegligible values are induced by 
renormalization group (RG) running. The approximate solutions to the RG equation for 
the left squark mass squared becomes 

imLL)ijif^ = M^eak) ^ {VcKM)^i (V'cKM)3j (^"^O + «o) log(- ^* 



.„2 * v''iJi^My3j v''CKMy3j ^■^""0 I "-o; "^sw^ )■, 

-'^-'weak 

and one can estimate {5'li)23 — 9 x 10"^. Note that this parameter is real however, since 
V23 and V33 are real in the Wolfenstein parameterization. Also other mass insertions are 
further suppressed by additional factors of and Y]^, so that I^^lrrlI^zI ~ 10~^ and 
1(5^^)23] ~ 10~^. These are too small to explain a large shift of 5',^;^ relative to S^k- 

In some SUSY GUT scenarios, the large mixing in the neutrino sector induces a large 
mixing in the right squark sector through the analogous RG running. The resulting mixing 
is about l2f 



Therefore the RR mixing is generically larger than LL mixing in SUSY GUTs where the 
see-saw mechanism is generating the nonvanishing neutrino masses. Still it cannot affect 
S(j}K by a large amount as discussed in the previous section, although the size of (5^ij)23 
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could be ~ 0(0.5) with a new complex phase so that B^-B^ mixing could be significantly 
shifted. 

On the other hand, if the SUSY fiavor problem is resolved by the alignment mecha- 
nism using some spontaneously broken fiavor symmetries, or decoupling (the effective SUSY 
models), the resulting LL or RR mixings in the 23 sector could be easily order of ~ 
(A = sm9c ^ 0.22) with either LL » RR, LL « RR or LL ^ RR ^ El], whereas the 
mixing in the 13 sector is further suppressed by additional power(s) of A('s) in order to avoid 
large contributions to B^-B^ mixing. These parameters will carry CP violating phases in 
general and can contribute significantly to AMg and sin2j3s- 

Further, in the presence of an LL or RR insertion, there can be an induced LR or RL 
insertion when_utan/3 is large with respect to the scalar masses, due to the double mass 



insertions 



Mm 



{hR)23 = (^Ll)23 X ^ . 



In case {^ll rr)'^'^ ~ 10^^, one can achieve {^f^RRLYit ~ 10"^, if /itan/? ~ 30TeV. This 
could be natural if tan/3 is large ~ 40. For larger LL,RR mixing, even smaller /itan/3 
would suffice to induce the needed LR, RL mass insertions of a size 10~^ — 10"'^. Since the 
Sll,rr.^s in SUSY fiavor models are generically complex, the induced (^i/^)23^ could carry 
a new CP violating phase js^], which can again explain a large negative Si^k- Also, the 
RR insertion can have a phase inherited from the neutrino mixing matrix in a SUSY GUT 
model, which would be transferred to the induced RL insertion j26|. Therefore the induced 
RL insertion can explain the observed S^k in SUSY GUT models. On the other hand, the 
phase of {S'[j^)23 for the case of minimal SUGRA is given by the CKM matrix elements, 
and it is real for the 2-3 mixing. Therefore the induced (5^^)23^ in SUGRA-like models is 
incapable of explaining the negative S^px- 



IX. CONCLUSIONS 

Recent data has given us hope that there may be new physics lurking in rare 5-decays, 
particularly B — > (pKs- In this paper, we considered several potentially important SUSY 
contributions to this process in order to see if a significant deviation in its time- dependent 
CP asymmetry S^k could arise from SUSY effects. In particular, we considered the SUSY 
gluino contributions in models with non-minimal flavor violation, and the Higgs-mediated 
contributions in models with minimal flavor violation. While the latter has the advantage 
of only affected the b — > sss transition, current bound on Bs fJ'^fJ'~ at Fermilab constrain 
the relevant operators to be too small. Models based on the LL and RR insertions in gluino 
interaction vertices also give contributions too small to alter S^k very much, unless gluinos 
are very light (close to the experimental bounds). However, in an LL or RR scenario one 
expects an impending observation of gluinos and squarks, and absence of B^-B^ oscillations 
at the Tevatron, since AMg increases dramatically in the LL (RR) case when S^k < 0. 

On the other hand, gluino-mediated LR or RL contributions can generate sizable effects 
in StpK as long as |(^L/j.ijL)23| ~ 10~^ - 10~^. As a by-product, we found that nonleptonic B 
decays such as B (j)Ks are beginning to constrain \ {S'[rrl)23\ as strongly as -B ^ Xs'j. We 
also studied various correlations among S^px, C^k, the direct CP asymmetry in i? ^ ^s1, 
AMg, All and sin2/3s. Using the LR or RL insertions, it is easy to obtain a negative Sij,K 
without conflict with any other observables. Furthermore, there are deflnite correlations 
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among S'^x, C(f,K and ^cp*'*^, and our explanation for the negative S^k can be easily tested 
by measuring these other correlated observables. In particular C^k can be positive only for 
the RL insertion. For this case, the direct CP asymmetry in 5 — > vanishes (assuming 
no RR insertion is present), unlike the case of the LR or LL insertions. In a scenario with 
both RL and RR insertions, the resulting direct CP asymmetry in i? ^ Xg'^ could be as 
large as in the LR mixing case, and there could be a large, complex contributions to 5°- 
5° mixing, leading to significant effects in AM^ and sin2/3s. The effects of the (pure) LR 
and RL insertions on B^-B^ are rather small and it would be difficult to distinguish our 
model from the SM by the B^-B'^ mixing, considering various theoretical and experimental 
uncertainties. 

Thus, we have found that three classes of supersymmetric models could not naturally 
explain a negative time dependent CP asymmetry S^px in i? — > (pKs if the initial hints 
of such an asymmetry should persist when data improves: LL, RR insertions, and Higgs- 
dominated decays. Two other classes, LR and RL insertions, can explain the data and 
can be distinguished to some extent with more and better data. Insertions of the size and 
kind needed can be generated naturally in simple string-motivated models and SUSY flavor 
models. 
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